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Abstract 

The search for the Higgs boson at the Tevatron and the LHC relies on detailed calculations of 
the kinematics of Higgs boson production and decay. In this paper, we improve the calculation of 
the distribution in transverse momentum, Qt, of the Higgs boson in the gluon fusion production 
process, gg — ► H, by matching the resummed distribution at small Qt with the C(a 4 ) fixed-order 
perturbative calculation at high Qt in the ResBos Monte Carlo program. The distribution is higher 
at large Qt than with the old 0(a 3 ) fixed-order perturbative calculation, and the matching with 
the resummed calculation is much smoother. The total cross section is also increased, more in line 
with next-to-next-to- leading-order calculations. We also study the effect of the new calculation on 
the distribution of Acftu in the overall process gg — > H — > W + W~ — > £ + £~ui', and the effect of 
PDF uncertainties on the distributions at the Tevatron and the LHC. 



1 



I. INTRODUCTION 



It has been over 25 years since the discovery of the W and Z gauge bosons in the UAl 
and UA2 colliders at CERN, and we are finally on the verge of discovering the source 
of their mass. With the on-going studies at the Tevatron collider at Fermilab and the 
turn-on of the Large Hadron Collider (LHC) at CERN, we will finally be able to probe 
directly the physics that breaks the electroweak symmetry and distinguishes the massive 
W and Z bosons from the massless photon. The simplest model of Electroweak Symmetry 
Breaking (EWSB) is the Standard Model (SM) of particle physics, which contains a complex 
electroweak scalar doublet that acquires a vacuum expectation value, thereby breaking the 
electroweak symmetry spontaneously. Three of the degrees of freedom of this complex 
doublet become the longitudinal modes of the massive W + , W~, and Z bosons, while the 
remaining degree of freedom is manifested as a single neutral scalar — the SM Higgs boson 
(H). Although this is not the only possible mechanism for EWSB, the search for the Higgs 
boson is the benchmark study for EWSB physics to be undertaken at the Tevatron and the 
LHC. If all goes well, it should be observed or ruled out within the next few years. 

The most stringent limits on the Higgs boson mass, m#, come from direct searches for 
the particle at LEP2 in the process e + e~ — > ZH, where the lower bound of 114.1 GeV has 
been obtained at 95% confidence level (C.L.) In addition, preliminary results from a 
combined fit of CDF and D0 data at the Tevatron has been used to exclude the mass range 
of 160 GeV < hih < 170 GeV at 95% C.L. [2] Beyond the direct search for a real Higgs 
boson, the effect of virtual Higgs bosons in loop calculations can be used to obtain indirect 
bounds on w,h- Current global fits to electroweak precision measurements, in combination 
with the direct search limit, prefer run < 191 GeV at 95% confidence level [3|. The Tevatron 
collider has a reasonable chance of discovery or exclusion over much of this preferred range of 
rax, assuming 7 fb _1 of data obtained by the end of Tevatron running [4|. Furthermore, the 
LHC can be expected to cover the entire range of Higgs boson masses up to about 1 TeV, 
which is a rough upper bound on m#, based on triviality and unitarity of the Standard 
Model fl. 

At both the Tevatron and the LHC the largest channel for production of the Higgs boson is 
gluon fusion, with the ggH coupling arising via a top quark loop. Other important channels 
are production of the Higgs boson with an associated W boson, Z boson, or top quark pairs, 
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as well as production of the Higgs boson through vector-boson or bottom-quark fusion. The 
importance of each channel for discovery /exclusion of the Higgs boson depends on its mass. 
A light Higgs boson (m# ;$ 135 GeV) decays predominantly to bottom quark pairs. In this 
case the inclusive Higgs boson signal is very difficult to pick out from the very large QCD 
bb background. The best sensitivity at the Tevatron in this mass range is instead found in 
the WH and ZH associated production channels, where the extra particles can be used to 
better distinguish the signal from background. A heavier Higgs boson (m# > 135 GeV) 
decays predominantly to W boson pairs with one of the W's potentially off-shell. In this 
mass range the inclusive production through gluon fusion is most important at the Tevatron, 
with the best sensitivity occurring around m# m 160 — 170 GeV, where the WW decay mode 
is fully open. These modes and many others were used in the recent combined fit at the 
Tevatron to exclude the mass range of 160 GeV < mg < 170 GeV at 95% C.L. |2|] Other 
important Higgs decay modes, both at the Tevatron and the LHC, are ZZ for high mass 
Higgs bosons, and tau pairs and photon pairs for low mass Higgs bosons. 

In order to best discern the Higgs boson signal from background, it is necessary to have 
the most accurate predictions possible for the kinematic distributions of the Higgs boson. 
In the leading order (LO) calculation of the gg — > H + X cross section, the Higgs boson is 
produced with exactly zero transverse momentum, Qt = 0. In higher order calculations it 
can have non-zero Qt, due to the emission of additional gluons or quarks, but the calculation 
at any fixed order of perturbation theory diverges as Qt — > 0. Thus, any fixed-order 
calculation is unsuitable for the study of the Qr-dependence of the Higgs boson (except 
at large Qt), or for the study of any other kinematic distribution that is strongly affected 
by soft gluon radiation. Fortunately, the soft-gluon effects that occur for small Qt can 
be incorporated into the calculation, either by their direct production in a parton shower 
Monte Carlo, such as PHYTHIA 6J or HERWIG Ml, or by analytic resummation of the 

nnn 

associated large logarithms, as proposed by Collins, Soper, and Sterman (CSS) |8|, |9j, [10]. 



This systematic resummation in powers of the strong coupling a s times powers of the large 
logarithm Iti(Q/Qt) has been applied to the Higgs boson process, as well as other processes, 



in the general resummation code ResBos 



111 ]. For the present process, the scale of the 



resummation, Q, is equal to the invariant mass of the produced Higgs boson, unless otherwise 
specified. 

The resummation of large logarithms at small Qt has been analyzed for Higgs production 
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in a number of studies in recent years 
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3, 
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13, uM- The power of the logarithms 



that are resummed is determined by parameters, which can be extracted order-by-order in 
a s from the perturbative Higgs production cross sections. The calculation of the Higgs bo- 

eading order, 
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son production cross section in gluon-gluon scattering has been calculated at 
next-to-leading order (NLO) [l9|] and next-to-next-to- leading order (NNLO) 
the infinite-top-quark-mass limit, and at LO and NLO 23j, |24| with full top quark mass de- 
pendence. In addition to the QCD corrections, the NLO electroweak (EW) corrections have 
also been considered in the infinite-top-quark- mass limit 25), and more complete calcula- 
tions have been performed by including light quark and top quark effects |26l. 1271] . Recently, 



the effects of the combined QCD and EW corrections were analyzed 



28, 291. The inclusive 



differential cross section at non-zero Qt, which begins at one higher power of a s , has been 



calculated at NLO jn the infinite-top-mass limit 
mass dependence 
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30. 
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321 ]. and at LO with full top quark 



34j. In the infinite-top-quark-mass limit, the heavy top quark loop 



contracts to an effective gluon-Higgs operator, which simplifies the calculation greatly, ef- 
fectively reducing the number of loops by one. In addition, it has been shown, at least at 
NLO, that it is a good approximation to calculate in the infinite top quark mass effective 
theory, while rescaling by the LO cross section with full top and bottom quark mass depen- 
dence {35)]. Thus, it has become standard to use this approximation to compute the Higgs 
boson cross section. For nonzero Qt this approximation is also good as long as Qt < m t 
and ran < rat [33J, [34 ] . 

A recent analysis using the resummation program ResBos to study the phenomenology 
of Higgs boson production at the Tevatron and the LHC was presented in Ref. JL7J. In that 
work the CSS resummation at small Qt was matched onto the LO calculation at large Qt- 
In this work we have updated the program so that it matches on to the NLO calculation at 
large Qt-, using the code developed in Ref. [32j. We shall see that this is more consistent 
with the precision currently included at small Qt in the resummation program. 

The remainder of the paper is organized as follows. In section [TT] we give a brief description 
of the resummation procedure that is implemented in ResBos. We explain the order of the 
coefficients used in the resummation calculation and how they are matched on to the fixed- 
order calculation at large Qt, and how the improvement is performed to include 0{ofy. 
We also calculate the total production cross sections of the Higgs boson at the Tevatron, 
the LHC with 10 TeV and 14 TeV center-of-mass (cm.) energy and compare with NLO 
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results. In section [TTT| we use the updated code of ResBos to produce various kinematic 
distributions for the Higgs boson. In particular, we emphasize the changes coming from the 
updated calculation, and note the improvement in the matching between the low and high 
Qt regimes. We also discuss the implications of the new predictions on the A0« correlation 
of the two charged leptons in the decay H — > W + W~ — > l + vi~v. Finally, in section HVl we 
give our conclusions. 



II. SOFT GLUON RESUMMATION 



A. Formalism 



In order to make transparent what we have implemented in the ResBos program, we 
begin by briefly reviewing the CSS formalism of soft gluon resummation. Resumming the 
soft gluons and using the narrow width of the Higgs boson to factorize the Higgs production 
from its subsequent decay, we can write the inclusive differential cross section for gg — > H — > 
VV — > 4 fermions as 

da(gg - HX - VVX - hhhhX) 



na 



Q 2 l 



M{H^VV^hf 2 hU) 



x 



(2 



2S {Q 2 -miy + (QiT H /m H y 

i- j d 2 b e lQT - b W gg {K, Q, x u x 2 , C 1X3 )W g N g p (b, Q, x 1 ,x 2 ) + Y(Q T , Q, x u x 2 , C 4 



(1) 



where S is the square of the center-of-mass energy; V and fi denote vector boson and 
fermion, respectively; Q, Qt, y, <Ph and T# are the invariant mass, transverse momentum, 
rapidity, azimuthal angle and total decay width of the Higgs boson, respectively, defined 
in the lab frame; and dU^ represents the four-body phase space of the Higgs boson decay, 



defined in the Collin-Soper frame 



361 ] . In Eq. ([!]), the quantity 



00 



V2G 



(2) 
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arises as an overall factor in the infinite-top-quark- mass limit. We have multiplied this by 
an additional factor 

a LO (m t ,m b ,m c ) 

'' : (7 L 0(00,0,0) ' 



which takes into account the masses of the top, bottom, and charm quarks at LO. It has been 
shown that multiplying the NLO Higgs cross section in the infinite-top-quark-mass limit by 
the factor k is a good approximation to the full mass- dependent NLO cross section over a 
wide range of Higgs boson masses [351]. In Eq. (pQ), • -)\ 2 denotes the matrix element 



squared of the Higgs boson decay whose analytical expressions are given in Ref. [17]. 

In Eq. (CQ), the term containing W gg dominates at small Qt, growing as Q^ 2 times a 
resummation in powers of In Q 2 jQ\, to all orders in a s . It can be expressed as 

W gg (b,Q,x u x 2l C 1 ^) = e- s ( b ' Q ' c ^J2( C 9"® f°)M(C gb ® h){x 2 ) , (4) 

a.b 

where the Sudakov exponent is given by 



S(b,Q,Cx,C 2 ) 



~~2 



Cl/b 2 P 



C\Q 2 



A(a.(fi), Cx) In ( ) + B(a.(fi), C x , C 2 ) 



(5) 
(6) 



The coefficients A and B and the functions C ga can be expanded as a power series in a 



and 



A(a s (fl),Cx) = £ (^-) n A^(Cx), (7) 

n=l ^ 71 ' 

B(a s (fi),Cx,C 2 ) = jr f^MY B^(Cx,C 2 ), (8) 

71=1 ^ 71 ' 

C ga (x,b,»,Cx,C 2 ) = £(^) n c£)(z,6,/i,^), (9) 

n=0 ^ ^ 

with fi = Cs/b. These quantities can be extracted order-by-order from the fixed-order 
calculations. In our numerical results, we have included A^ l,2 ^\ B^ 1 ' 2 ^ and C^ Q ' l \ whose 
analytical expressions are given in Appendix |A] for completeness. We use the canonical 
choice for the renormalization constants, C x — C% — 2e~ 7B , C 2 = C± = 1, which simplifies 
the above expressions. The function W?L P describes the non-perturbative part of the soft- 



gluon resummation, in which we use the BLNY parameterization [37]. 

Finally, in Eq. (jTJ), the term containing Y incorporates the remainder of the cross section, 
which is less singular as Qt — > than the H^-term. It consists of the difference between the 
full cross section at finite Qt and the small-Qr limit of this cross section, each calculated 
to the same order in a s . At small Qt, these cancel, so that the contribution of the F-term 



is small, and the resummed term dominates. At large Qt, where the logarithms become 
small, the resummed term cancels against the small-Qr limit term (to the given order in 
a s ), so that the cross section approaches the fixed-order calculation. More details of how 
this matching process between the resummed calculation and the fixed-order calculation 
is implemented in ResBos can be found in Ref. In previous studies of Higgs boson 

production at hadron colliders [ijl, 1^], the high Qt perturbative calculation was included 
in ResBos at 0(a^). The major update to the program that we have incorporated in this 
paper is to include the high Q T perturbative calculation at This was done by using 

the code of Ref. |32J to rescale the perturbative piece of the grids used by the ResBos code by 
the factor (Pert(a^)+Pert(a^))/Pert(a^), where Pert(a^) and Pert(o^) are the contributions 
to the Qt distribution of the Higgs boson of order and a^, respectively. 

Note that the description of the perturbative piece in terms of LO, NLO, or NNLO is 
problematic; for example, 0{a\) would be considered as NLO when referring to the (non- 
zero) transverse momentum distribution, but it would be considered as NNLO when referring 
to the total Higgs production cross section. Thus, we will refer to the power of a s when 
comparing the precision of the perturbative piece of the calculation used in this work (0(o4)) 
versus that used in the previous works (0(a^)). 



B. Predictions of the total cross section for gg —* HX at hardon colliders 

The primary use of the resummation code ResBos is for the calculation of the transverse 
momentum spectrum, as well as other distributions that are influenced strongly by soft gluon 
effects. However, it also gives a calculation of the total cross section that is comparable 
to that of a fixed-order calculation, depending on the order to which the resummation 
coefficients have been included. We have included all of the coefficients {A^ ,2 \ C^ ' 1 )) 
in ResBos that are necessary to produce a NLO calculation of the cross section. In addition, 
we also have included the NNLO coefficients and and are only missing the function 
that is necessary to give a NNLO calculation of the total cross section. The function 
should be extractable from the NNLO analytic expression of the cross section, but this 
has not been achieved as yet. Thus, our calculation of the cross section should be comparable 
to a NLO fixed-order calculation, and in fact contains much of the (presumably dominant) 
contributions at NNLO. 
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LHC 14 TeV 


LHC 10 TeV 


Tevatron 1.96 TeV 


m H (GeV) 


RES 


NLO^ 5 


NLO m g iu 


RES 


NLO RES 


NL0 Higlu 


RES 


NLO RES 


NLQ Higlu 


140 


32.12 


25.79 


27.23 


18.09 


14.37 


15.16 


0.50 


0.37 


0.39 


150 


28.85 


22.97 


24.09 


15.85 


12.67 


13.28 


0.40 


0.30 


0.32 


160 


25.61 


20.61 


21.50 


14.05 


11.25 


11.74 


0.33 


0.24 


0.26 


170 


23.03 


18.60 


19.31 


12.54 


10.06 


10.45 


0.27 


0.20 


0.21 



TABLE I: Production cross section of the SM Higgs boson via the gluon fusion process, gg — > HX, 
in pb at the LHC and Tevatron. We show the results for both 14 TeV and 10 TeV cm. energy 
at the LHC and 1.96 TeV cm. energy at the Tevatron. The label RES indicates the results of 
resummation calculations predicted from ResBos, and the labels NLO^ 5 and NLO^* 9 '" indicate 
the NLO results calculated from ResBos and Higlu codes, respectively. 



In the remainder of this section we compare the predictions for the total cross section 
from ResBos against NLO predictions. In table [H we present the total cross sections of 
gg — > HX for several benchmark points from the resummation (RES) calculations using the 



updated ResBos program with CTEQ6.6M PDF 38]. These are compared to an expansion 
of the resummation formula to NLO (RES-NLO), and also to an exact NLO calculation. The 
former is an exact NLO QCD calculation with the same implementation for including the 
effect of the masses of the top, bottom, and charm quarks at LO, cf. Eq. (3). The latter is 
calculated with the help of the public code HIGLU 39J.* 1 We note that the main difference 
between these two NLO calculations is in the handling of the quark mass dependences. As 
explained in section III Al the resummed calculation, as well as the RES-NLO calculation, 
is performed in the heavy top quark mass limit, with the quark mass dependence included 
by the LO factor k given in Eq. [3j On the contrary, HIGLU uses the exact NLO two-loop 
calculation, including both the top quark and the bottom quark in the loop. 

The cross section is consistently higher in the resummed calculation (RES) than for the 
NLO calculations, due to the enhancement from the NNLO corrections. This is seen more 
easily in Fig. fTJ which displays the cross section of the Higgs boson production as a function 
of tuh at the LHC and Tevatron, where the black, blue and red curves denote the RES, 



# 1 More detailed analysis of the HIGLU calculation, e.g. the PDF uncertainties and scale dependence, are 
given in the Appendix iBl 
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FIG. 1: The resummation the total cross section for gg — * HX at the LHC and the Tevatron 
using the updated ResBos program with CTEQ6.6M PDF. We consider both 10 TeV and 14 TeV 
cm. energy at the LHC and 1.96 TeV cm. energy at Tevatron. The NLO predictions from ResBos 
and Higlu programs are also shown. 

RES-NLO and HIGLU results, respectively. We display explicitly the enhancement in the 
cross section in the resummed calculation in Fig. [2j where we plot the ratio of the RES 
cross section to the NLO cross section, both for the RES-NLO and the HIGLU calculations. 
The ratio of RES to NLO drops rapidly when van ^ 300 GeV and reaches a minimum for 
mil ~ 380 GeV. This unusual dependence on m# can be traced to the handling of the 
virtuality of the Higgs boson in the calculations. The HIGLU calculation is for an explicitly 
on-shell Higgs boson production without the subsequent Higgs boson decay, i.e. 5{Q 2 — m 2 H ). 
The same is also true for the RES-NLO calculation. On the contrary, the ResBos code takes 
into account the Breit- Wagner width effects, see Eq. (CQ), to perform a realistic simulation. 
Thus, the cross section is calculated for an off-shell Higgs boson of mass-squared Q 2 , which 
is then convoluted with the Breit- Wigner. This gives a sizable effect when the Higgs width is 
large, and in particular for masses around 300-400 GeV {m H > 2m t ) where the cross section 
curves have noticeable structure. 

III. NUMERICAL STUDY OF gg -> H -» WWW -> ttvv 

For our analysis of kinematic distributions, we will use several benchmark values for the 
Higgs boson mass: m^j = 140, 150, 160 and 170 GeV at the Tevatron RUN II, and ran = 160 
GeV at the LHC. We will focus on the most promising discovery mode for this mass range, 
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m H (GeV) m H (GeV) 

FIG. 2: The ratio of the total cross section of gg — * HX predicted from the resummation calculation 
to that from NLO calculations at the LHC and the Tevatron. 

which is gg —>■ H — > W + W~ — > £ + £~b>v, where denotes a charged lepton and v{y) is a 
neutrino (anti- neutrino). Note that all of the cross sections given in this section include the 
decay branching ratio of the Higgs boson and that the flavors of leptons are not summed 
over. 

We begin with the transverse momentum distribution of the Higgs boson at the Tevatron, 
displayed in Fig. [31 The cyan lines and the black dash lines are the predictions with matching 
at high Q T to 0{af) and 0(aP s ) fixed-order perturabative calculations, respectively, and the 
bands show the PDF uncertainties. The peak position (at Qt ~ 10 GeV) is the same 
for both the dash and cyan lines, since this is determined exclusively by the resummed 
contribution to the calculation. At high Qt-, however, the distribution is determined mainly 
by the perturbative contribution, and it is substantially higher for the 0{af) calculation. 
We also note that the 0(a\) perturbative contribution matches much more smoothly with 
the resummed contribution. This is particularly apparent when plotted with a logarithmic 
scale, as in the lower panel of Fig. [3j where the unphysical kink around Qt ~ 70 GeV in the 
0{pt Z g) curve is absent in the 0(oq) curve. Similarly, we show the improved Qt distribution 
of the Higgs boson at the LHC in Fig. HI including PDF uncertainties as well. The kink 
seen in the 0(a^) dashed line is shifted to higher Qt (~ 110 GeV) compared with the case 
at the Tevatron. With matching to the 0{oQ) calculation, the curve becomes much more 
smooth. The uncertainty at the peak region is about 8% and increases to higher than 10% 
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FIG. 3: Transverse momentum, Qt, of the SM Higgs boson at the Tevatron. The plots in the 
lower panel are the same as those in the upper panel, except for using a logarithmic scale for the 
y-axis. The dashed lines and cyan lines are the calculations when matching at high Qt to the 
0(a^) and O(af) fixed-order perturbative contributions, respectively, and the bands reflect the 
PDF uncertainties.. 

when Qt is larger than 40 GeV for all curves at the Tevatron, while the uncertainty stays 
at about 2% ~ 6% for the full range of Qt at the LHC. 

We can use the updated Qt distributions at the LHC and the Tevatron to study the 
dependence of this distribution on the collider cm. energy. The peak position, Q^(peak), 
only changes mildly with the cm. energy. For example, Q^peak) is roughly 10 GeV, 13 
GeV and 13 GeV at the Tevatron, LHC 10 TeV and LHC 14 TeV cm. energy, respectively, 
for run = 160 GeV. On the other hand, the average transverse momentum, (Qt), increases 
more substantially when the mass of the Higgs boson or the energy of the collider increases. 
For rrtu — 160 GeV and taking the average in the region of < Qt < 200 GeV, we find 
that (Qt) increases from about 26 GeV at the Tevatron to about 40 GeV and 43 GeV at 
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FIG. 4: The same as Fig. d but for the SM Higgs boson at the LHC. 
70 i 1 1 1 1 1 1 1 1 1 1 40 i — i — i — i — i — i — i — i — i — i- 



J i I i I i L 



Tevatron 1 .96 TeV 




_i i i i_ 



100 200 300 400 500 600 100 

m H ( GeV) 



150 

m 



200 



H (GeV) 



FIG. 5: The average of Q T (in the region of < Q T < 200 GeV) of the Higgs boson at the LHC 
and the Tevatron. 

the LHC with 10 TeV and 14 TeV cm. energy, respectively. The dependence of (Qt) on 
the Higgs boson mass at the various colliders is shown in Fig. [5j 

Finally, we consider the distribution in the difference in azimuthal angle of the two charged 
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leptons in the Higgs boson decay, A<j>u, which is useful for making cuts to extract the Higgs 
boson signal from background. This distribution is shown for the Tevatron in Fig. [6j The 
solid and dash lines in the lower panel are the calculations when matching to 0{oq) and 
0(al) at high Qt, respectively. From these plots we see that the signal process gg —>■ H — ► 
W + W~ — > i + i~vv prefers to have the charged leptons both moving in the same direction 
near A0« « 0. This can be understood from angular momentum conservation and the 
left-handed nature of the PF-boson decays. In the Higgs boson rest frame, the W + and W~ 
are produced back-to-back with opposite spins. The W~ boson decays with the charged 
lepton momentum anti-correlated with the W~ spin, while the W + boson decays with the 
charged anti-lepton momentum correlated with the W + spin. As a result the two charged 
leptons tend to move in the same direction. This feature still holds even when only one of 
the W bosons from the Higgs boson decay is on-shell. On the other hand, the background 
events, which predominantly originate from W + W~ pair production, are more likely to be 
produced with the charged leptons back-to-back with A</>« « tt. 

The enhancement of the distribution due to the improved calculation is exhibited in the 
plots in the upper panel of Fig. [6J where the blue line is the ratio of the solid curve to the 
dash curve. From these figures we see that the new calculation with matching at 0(aj) 
enhances A0« by about from 14% to 18% over the old calculation. Interestingly, the largest 
enhancement occurs at A0« « 0, while the enhancement has a minimum around A0« « 2 
radians. Similarly, in the A0« distribution at the LHC, which is shown in Fig. [7J the 
enhancement due to the improved matching varies between 16% ~ 23%, which is slightly 
larger than for the Tevatron. The change in the shape of the A0« distribution between 
the 0(a%) and the 0(a\) calculations can be considered purely kinematical, since the Higgs 
boson is spinless, and the comparable distribution in the Higgs boson rest frame would 
not be affected by changes in the production cross section. In the lab frame, however, the 
distribution is affected by the fact that the cross section is larger at high Higgs boson Qt in 
the 0(aig) calculation, so that more of the charged lepton pairs are produced with a bigger 
momentum boost. 



13 



I I I I I I I I I I I I I I 

140 GeV - 




i i i i I i i i i I i i i i 



12 3 
A<|) ;( (rad) 



.1 i i i I i i i i I i i i i 
140 GeV 



i i i i I i i i i I i i i i 



R 

1.18 
1.17 
1.16 
1.15 
1.14 

3 



i i i i I i i i i I i i i i I 
150 GeV - 



Tevatron RUN II 
R 




i i i i | i i i i | i i i i |_ 
160 GeV - 



1 2 3"'"0 
A§ n ( rad) 



R 

1.18 
1.17 
1.16 
1.15 
1.14 
1.13 



.i i i i I i i i i I i i i i 
170 GeV " 



i i i i I i i i i I i i i i 



12 3 

A§ u (rad) 







. i i i i | i i i i | i i i i |. 
150 GeV 



I i i i i I 




1 2 3 
A(|) /( f rad) 







. i i i i | i i i i | i i i i |. 
160 GeV 



I i . i i I 



1 2 3 
A(|) ;; ( rad) 



.i i i i | i i i i | i i i i 
170 GeV 



I i i i . I 



1 2 3 1 2 
A(|) /( ( rad) A(|) ;/ ( rad) 



1 2 3 
A(f> ;; ( rad) 



FIG. 6: Acf)££ distributions at the Tevatron. The solid and the dash curves in the lower panel 
are the calculations when matching at high Qt to the 0(a^) and 0{ot^) fixed-order calculations , 
respective rve. 
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FIG. 7: The same as Fig. [6] but for the LHC. 
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IV. CONCLUSION 



The search for the Higgs boson at the Tevatron and the LHC relies on detailed knowledge 
of the distributions (event rate and shape) of the Higgs decay products. These, in turn, are 
sensitive to the kinematics of the Higgs boson production. In this paper, we improve the 
calculation of the transverse momentum, Qt, distribution of the Higgs boson in the gluon 
fusion production process, gg — > H, at the Tevatron and the LHC by matching the resummed 
distribution at small Qt with the 0(aj) fixed-order perturbative calculation at high Qt in 
the ResBos Monte Carlo program. The total cross section of gg — > HX predicted from the 
updated ResBos is always larger than NLO calculation and the enhancement can reach about 
40% at the LHC for a heavy Higgs boson. The difference between RES and NLO calculation 
becomes minimum (~ 10%) in the threshold regime of > 2m(. The unphysical kink, 
which had been obtained when matching at high Qt with the 0(a z s ) fixed-order perturbative 
calculation, is removed. The PDF uncertainties are studied in the Qt distributions as well. 
At the Tevatron we find the uncertainty to be about 5% in the peak area and larger than 
10% in the region of Q T > 50 GeV. At the LHC, the uncertainty is about 2% ~ 6% for 
the entire region of Qt- For the comparison of Qt, we found that the average value of Qt 
increases when the cm. energy or the mass of the Higgs boson increases, however, peak 
position is insensitive to both. We finally study the opening azimuthal angle , A0«, of the 
two charged leptons in the Higgs boson decay, H — > W + W~ — > l + l~vv. The distribution of 
A0« is enhanced by 14% ~ 18% and 16% ~ 23% at the Tevatron and the LHC, respectively, 
with the largest enhancement occurring at small opening angles. With current and expected 
integrated luminosity, the Tevatron has the capability to discover or exclude the SM Higgs 
boson over a significant range of mass. Detailed predictions of the kinematics of the Higgs 
boson, such as its transverse momentum distribution, play a crucial role in this analysis. 
This is our motivation for the continued improvement of the ResBos program and for this 
study. 
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APPENDIX A: A, B, AND C COEFFICIENTS 



For completeness, we give expressions for the so-called A, B, and C coefficients used 
in our resummation calculations. In our numerical result, we have used A^ 1,2 \ B^ 
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42|], |43| and 



.i) 



44 



4Q|, 



45| . Their analytical expressions for the process 



gg — > H are much simplified in the canonical choice of the renormalization constants, C\ = 
C3 = 2e _7E , (72 = 6*4 = 1, which we use in this project. For this choice of renormalization 
constants, we have: 
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where (3 = (11N C - 2N f )/Q and T R = 1/2; 

C$f{*) = S(l-x), (A6) 

Cj»(x) = 0, (A7) 

Cj£(x) = 6(1 -x)^-^, (A8) 



CftXx) = ?fx, (A9) 



where x is the momentum fraction carried by the gluon after splitting from its mother 
particle (gluon g or quark q). 

APPENDIX B: TOTAL CROSS SECTIONS OF gg -> HX AT NLO 

In this section we evaluate the uncertainties in the cross section due to uncertainties 



in the PDFs and due to higher order corrections, as i 
dependence. Since we have used the code HIGLU 



luminated by renormalization scale 
as a reference comparison for the 
total cross section in our resummation calculation (see for example Fig. [1]), we will use it 
to produce numerical results in this section. We expect the uncertainties in the total cross 
section in the resummed calculation to be comparable. Note that the Higgs boson decay 
is not implemented in HIGLU, so all the results presented in this appendix are for on-shell 
Higgs boson production only. 

In Fig. [HI we display the uncertainties both due to the PDF uncertainties and due to the 
scale dependence. The uncertainties due to the PDFs (primarily the gluon PDF), relative 



to the cross section with the best-fit PDF, is shown by the bands in 



'ig. [HJ The PDF error 

is calculated from the master formula given in Eq. (2.5) in Ref. 46j, using the 44 sets of 
CTEQ6.6M package. For 100 GeV < hih < 600 GeV, the uncertainty is smaller than 5% at 
the LHC with a cm. energy of both 14 TeV and 10 TeV. In the intermediate mass region, 
200 GeV < iiih < 300 GeV, where the gluon PDF is more constrained, the uncertainty is 
reduced to about 2% ~ 3% Setting x\ ~ = x where x\^ is the momentum fraction of the 
incoming parton, and using s = XiX 2 s, we obtain (x) ~ rn H /^/~s. From this we see that the 
minimum in the PDF uncertainty in both Fig. [8(a) and (b) occurs around (x) ~ 0.022. On 
the contrary, the Higgs boson production cross section at the Tevatron suffers from much 
larger PDF uncertainties. For example, the uncertainty increases from 5% to 14% for the 
mass of Higgs bosom mass 100 GeV < mn < 200 GeV. 
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FIG. 8: The PDF uncertainties of the NLO total cross section of gg —* HX shown in bands and the 
dependence on the renormalization scale at the LHC and Tevatron using the CTEQ6.6M PDFs. 
The renormalization and factorization scales (a = (ar = /j>f are set to /Uo/4, fio/2 and 2/xo, where 
A*o = m Hi an d the ratio is taken with respect to the cross section evaluated with fi = fiQ. 

We also display the uncertainties in the cross section calculation at NLO due to the renor- 
malization scale (fiR) and factorization scale (fip) dependence in Fig. El These uncertainties 
can be considered as an estimate of the size of the unknown higher order corrections. For 
this study, we have set fi = fiR = \if and vary it around the central value of /i = rn H . Typ- 
ically, a factor of 2 is used to estimate the size of the higher order corrections, so we have 
displayed curves with /i = 2/i n and n = jio/2. In addition, since the NNLO QCD corrections 
prefer a scale of fi = /Uo/4 [21], we also display a curve with that value. In Fig. E]we plot 
the ratio a(ni)/a(no) as a function of tuh both at the LHC and at the Tevatron. The cross 
sections vary between about —15% for /z = 2/i and +20% for /i = fi /2 at the LHC, and 
can reach about +40% when using /i = /io/4. At the Tevatron, the scale dependences are 
even larger. We note that the scale dependence at both colliders is insensitive to m# and 
dominates over the PDF uncertainties. 

The PDF uncertainties can be improved further by using the new set of CTEQ PDFs 



(named CT09 [47]), which take into account the recent inclusive jet data at the Tevatron [48, 
491 ] . Fig. [9] shows the relative error, 5a, in the Higgs boson production cross section due to 
PDF uncertainties derived from CT09 (black) and CTEQ6.6M (red). The uncertainties in 
the Higgs boson production cross section are improved substantially using the new set of 
PDFs, both at the LHC and at the Tevatron. Due to the modification of the gluon PDF, 
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FIG. 9: The uncertainties of the NLO total cross section of gg — ► HX due to the PDFs. The red 
curves are for CTEQ6.6M and the black curves are for CT09. 

the minima of the PDF uncertainties at the LHC are shifted to smaller values of the Higgs 
boson mass. 
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